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Abstract 

Based on observational tests and constraints on halo structure, dark matter is generally taken 
to be cold and essentially collisionless. On the other hand, given the large number of particles 
and forces in the visible world, a more complex dark sector could be a reasonable or even likely 
possibility. This hypothesis leads to testable consequences, perhaps portending the discovery of 
a rich hidden world neighboring our own. We consider a scenario that readily satisfies current 
bounds that we call Partially Interacting Dark Matter (PIDM). This scenario contains self- 
interacting dark matter, but it is not the dominant component. Even if PIDM contains only a 
fraction of the net dark matter density, comparable to the baryonic fraction, the subdominant 
component's interactions can lead to interesting and potentially observable consequences. Our 
primary focus will be the special case of Double-Disk Dark Matter (DDDM), in which self- 
interactions allow the dark matter to lose enough energy to lead to dynamics similar to those in 
the baryonic sector. We explore a simple model in which DDDM can cool efficiently and form a 
disk within galaxies, and we evaluate some of the possible observational signatures. The most 
prominent signal of such a scenario could be an enhanced indirect detection signature with a 
distinctive spatial distribution. Even though subdominant, the enhanced density at the center 
of the galaxy and possibly throughout the plane of the galaxy (depending on precise alignment) 
can lead to large boost factors, and could even explain a signature as large as the 130 GeV Fermi 
line. Such scenarios also predict additional dark radiation degrees of freedom that could soon 
be detectable and would influence the interpretation of future data, such as that from Planck 
and from the Gaia satellite. We consider this to be the first step toward exploring a rich array 
of new possibilities for dark matter dynamics. 

1 Introduction 

All known particles make up only a small fraction of the energy density in our universe, yet the 
Standard Model is extremely complicated: three forces, one Higgsed, one confining, plus quarks 
and leptons organized into three generations. This model — the components of the visible universe — 
deviates markedly from any apparent principle of minimality. Yet, when considering the 85% of 
the matter in the universe that is dark, our usual response is to turn to minimal models of a cold, 
collisionless particle: a WIMP, perhaps, or an axion. Slightly less minimal variations are sometimes 
studied, often motivated by data that is in mild conflict with the cold dark matter paradigm. Self- 
interacting dark matter (SIDM) [I] and warm dark matter [2] are well-studied examples. These 
minimal choices are, to some extent, justified by Ockham's razor. We know that the bulk of 
galaxy halos consists of dark matter organized into large, diffuse, spheroidal distributions, and 
(based on halo shapes and the Bullet Cluster) that the stuff making up these halos is approximately 
collisionless, which offers some support for the idea of minimality. Still, confronted with the richness 
of physics in the visible world around us, it is tantalizing to imagine that the dark world could be 
similarly complex, full of structures, forces, and matter that are invisible to us. We might hope that 



a whole sector of the universe as rich as our own exists just out of sight. Our goal in this paper 
is to argue that this superficially fanciful idea should be taken seriously as a testable hypothesis, 
which might even help to resolve some of the deficiencies of the ACDM scenario. Double-Disk Dark 
Matter (DDDM) is a concrete form of this idea, in which a small fraction of all dark matter has 
dissipative dynamics causing it to cool into a disk within the Milky Way galaxy. Double-Disk Dark 
Matter acts less like typical noninteracting dark matter than like a new kind of ordinary matter, 
constituting an invisible world that may be literally parallel to our own. 

The general scenario that we propose, Partially Interacting Dark Matter (PIDM), is that a sub- 
dominant component of dark matter has self-interactions. The more specific DDDM scenario that 
we focus on assumes the existence of a massless (or nearly massless) U(l) gauge boson that permits 
dissipative dynamics. This will generate observationally distinguishable dark matter consequences 
and in particular a thin dark matter disk similar to the baryonic disk when an additional light dark 
charged particle is present with sufficient abundance. Even though by assumption the dark matter 
is subdominant, the density of the interacting component in the disk could be much higher than the 
dominant diffuse dark matter that is spread throughout the halo and could hence lead to stronger 
indirect signatures. These could include observations such as the recently reported Fermi line (3j|4], 
which is challenging to explain without significant boost factors or tuned scenarios (many of which 
are summarized in [5][6]). We will construct models that generate this signal to illustrate one pos- 
sible consequence of DDDM, but the potential signatures span a much wider range of possibilities. 
Data will determine whether a richer dark matter sector does in fact exist. 

Because existing constraints are weak as we will see in Section [2j the fraction of additional 
interacting dark matter is fairly unconstrained so long as it is subdominant and nondissipative. 
When the dark matter does collapse into a disk, the constraint is stronger due to the Oort limit, 
which is a bound on the amount of matter in the vicinity of our solar system. For many of the 
numerical results in this paper, we will assume that the energy fraction is at the Oort limit, yielding 
a fraction of the energy density comparable to but smaller than that of baryons. Because new forms 
of matter might be present at the weak scale, and because our scenario was motivated in part by 
the Fermi photon line, we will often assume a new interacting dark matter particle of mass near 100 
GeV. Again, the scenario allows for a wide range of parameters and we consider this mass simply 
for specificity. 

Partially Interacting Dark Matter and the Double-Disk Dark Matter scenario that we focus on 
can potentially lead to other interesting consequences. In future work we would like to explore the 
effect of this component on overall structure formation and on observations such as gravitational 
lensing and detailed sky surveys. This scenario can also lead to an acoustic oscillation signature 
similar to that from baryons. Detailed observations of the CMBR and galaxy correlation functions 
should also detect or constrain DDDM. 

We note that other authors have proposed interesting scenarios for dark matter involving long 
range forces and bound states: the older idea of mirror matter pl]12|, recent work on dark matter 



charged under a hidden U(l) [13-15 , and more recent work on dark atoms |16f[l8] . In fact, all of 
these ideas were foreshadowed in much earlier work by Goldberg and Hall [19] . Furthermore, the 
idea of self-interacting dark matter (SIDM) has been studied intensively as a possible solution to 



discrepancies between cold dark matter models and observations (e.g. the cusp/core problem 20 



2l]), beginning with Ref. [T]|22j (for recent progress see Refs. |23[|24| ). These scenarios are often 
very constrained by observations of the halo and galaxy and galaxy cluster interactions. The 
chief difference in our scenario is that the interacting component of dark matter is subdominant. 
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We will see that such a scenario is far less restrictive since galaxy shape constraints as well as 
direct constraints on interactions rely predominantly on the existence of a dominant noninteracting 
component. A subdominant component can interact and permit much richer dissipative dynamics 
and that is what we consider below. 

We begin by explaining in Section [2] how current constraints can allow interesting amounts of 
interacting dark matter. In Section [3] we explain that they also allow background dark radiation in 
the amount predicted by our model, which will be probed in the near future by the Planck satellite. 
In section [4] we argue, using the tentative Fermi line as an example, that possible gamma ray line 
signals of dark matter annihilation at the level probed by current observations are difficult to explain 
without a very large boost factor. In section [5j we show that it is possible for a subdominant species 
of dark matter to efficiently cool within the lifetime of the universe. We follow this with a discussion 
of the structure that forms as a result of cooling, the DDDM disk, in section [6] At this point, having 
established that a disk can form and estimated how much matter can lie in it, we are equipped to 
return to the issue of a boost factor for indirect detection, which we discuss (along with direct 
detection prospects) in section [7] Even if it does not provide such signatures, DDDM could exist 
and be detected through its gravitational effects. We conclude in section [8] with an outline of the 
many exciting possibilities for future work on this subject. 



2 Constraining the Amount of Allowed DDDM Abundance 

We define the fraction of the energy density in PIDM, compared to ordinary dark matter, as 
e n = ^Ifp™ ■ Furthermore, we expect the relative fraction of different matter components in the 

Milky Way is comparable to that in the universe as a whole. So we take cq ~ ^ph^m/^dm wnere 
-^PIDM ^ s ^ e total mass of PIDM in the galaxy and is the total mass of all dark matter in 

the galaxy. In fact, in all likelihood the ratio can be bigger in terms of the total energy accounting, 
since our strongest bound is only on PIDM in disk form and only about a third of the baryons end 
up in the disk. For this case of DDDM, we denote the fraction of mass in the Milky Way's disk by 

~ M gal ' 

If DDDM is organized similarly to baryons, the total energy fraction in DDDM would be more like 
~ 3e. But in most of the paper we will take ~ e for simplicity. 

Current bounds on self-interacting dark matter arise from halo shapes and cluster interactions. 
So far such bounds have been calculated only for a single dark matter component, for which they 
can be quite constraining. Self-interactions lead to more spherical halos, especially in the inner 



region, where the density is higher and interactions are more frequent 22 23 . 

One bound of this type arises from the halo for the galaxy cluster MS 2137-23 25 , which is 
measured by gravitational lensing to show a 20% deviation from axial symmetry at radius 70 kpc 
from the center. Another bound comes from measurements of X-rays emitted by hot gas in the 
elliptical galaxy NGC 720 |26| , showing 35% deviations from sphericity at distances of 5 to 10 kpc 
from the center. Roughly speaking, these bounds on SIDM exclude the possibility that a typical 
dark matter particle has scattered at least once in the age of the universe. The time scale for a 
dark matter particle to scatter is (nav) -1 , with n = p/m the number density of other dark matter 
particles it could interact with. Hence, the limits are typically expressed as a bound on cross section 
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per unit mass, a/m. The inferred bound is approximately a/m < 0.1cm 2 /g ~ 0.2 barn/GeV 
Such cross sections are large by the standards of pointlike particles, but readily arise through long- 
range forces (for instance, Rutherford scattering has a l/v 4 enhancement at low velocities, leading 
to large cross sections) or through large composite objects (e.g., atoms with Bohr radii ao ~ m^ 1 /a 
at weak coupling a, which have cross sections even larger than AttQq 16 ). In the case that all 



dark matter consists of a particle (and its antiparticle) charged under a massless U(l) gauge boson, 
these bounds have been studied in Refs. |14) and [15] , which find that because of the low- velocity 
enhancement of Rutherford scattering the bounds exclude a thermal relic abundance. In the case 
of composite dark atoms, the bounds have been studied in Refs. 16 , 18 . They exclude a portion of 
the parameter space, but less than in the fully ionized case, because atom/atom scattering is closer 
to a hard-sphere interaction without a long-range force. 

Such bounds do not directly apply to PIDM since a sufficiently small fraction of all matter could 
have extremely strong interactions without affecting observations at all. A conservative estimate of 
the allowed abundance of PIDM can be found using the recent halo shape analysis of Ref. [24] , which 
corrects certain deficiencies in earlier analyses and argues that a/m = 1 cm 2 /g is ruled out by the 
X-ray observations of NGC 720 but that a/m = 0.1 cm 2 /g is allowed by all current bounds. Figure 
5 of Ref. [24] shows that 20% to 30% deviations from spherical symmetry, as observed in data, are 
compatible with all dark matter particles scattering a few times in the age of the universe. In that 
case, 10% of all the DM would have scattered 10 times or more. PIDM — which will scatter multiple 
times and form denser structures — together with 90% ordinary dark matter with no interactions at 
all, will leave the halo triaxial. Numerical simulations should be done to quantify these statements 
more carefully and produce a definite bound, but we summarize by saying that halo shapes are 
compatible with 10% of the dark matter having arbitrarily strong interactions if the other 90% 
does not interact at all. Even more may be allowed. One reason that this is plausible is that we 
already know that 15% of all matter interacts strongly and forms dense structures not captured by 
simulations of dark matter halos, without violating the constraints from observation: we refer, of 
course, to ordinary baryonic matter. 



Observations of the Bullet Cluster have also set bounds on dark matter self-interactions 27 28 



which for SIDM have been argued to be weaker than those from halo shapes 15 24] . They are, 
however, more readily interpreted as a bound on PIDM than the halo shape bounds. In the Bullet 
Cluster, two merging clusters have led to separation between collisional material (hot gas) and 
collisionless material (stars and ordinary dark matter). The mass of dark matter in the subcluster 
with stars (inferred from gravitational lensing) leads to the conclusion that no more than 30% of 
the dark matter has been lost to collisional effects. Thus we expect that, if the bulk of dark matter 
is completely collisionless, the Bullet Cluster bound tells us that a subdominant component making 
up 30% of all dark matter can have arbitrarily strong self-interactions. It would be left behind with 
the gas, without changing the lensing observation that tells us the dominant component of dark 
matter moved with the collisionless stars. 

In fact, the most stringent PIDM abundance constraint arises only when there is dissipation and 
a disk is formed. In that case, stellar velocities in and out of the plane of the galaxy yield stronger 
bounds. Such velocity distributions offer interesting prospects for reconstructing the galaxy's grav- 
itational potential and hence inferring the distribution of matter within it. The most relevant 
bounds to date come from the Oort limit, i.e. the inferred local density of matter near the Sun 
from observations of nearby stars. A recent determination that the local dark matter density is 



0.3 ± 0.1 GeV/cm 3 29 relied on the kinematics of stars between 1 and 4 kpc above the galactic 
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plane. Older results were based on stars within 100 pc of the Sun, surveyed by the Hipparcos 



satellite, and another sample of stars extending out to 1 kpc 30 , 31 1 . In the presence of a possible 



dark disk |32| these observations were estimated to be consistent with a local dark matter density 



between 0.2 and 0.7 GeV/cm 3 |3ll. For our purposes, the most convenient form of the bound is the 



constraint on the surface density measured below a height Zq, which is defined by: 

E(|;z| < z ) = / p(z)dz. (2) 



J-Zq 

£ is approximately equal to the vertical gradient of the gravitational potential, (2ttG d z <& , which 
determines the vertical acceleration of stars. As quoted in [31] , the total surface density inferred 
from stellar kinematics is £tot(M < 1-1 kpc) = 71 ± 6 M /pc 2 . The surface density inferred from 
visible baryonic matter (stars, stellar remnants, and interstellar gas) is S v ; s = 35 to 58 M /pc 2 . 
We interpret the difference between these numbers as an approximate measure of the amount of 
DDDM allowed by data. The ranges are one sigma error bars, from which we conclude that at 95% 
confidence level the amount of surface density in nonbaryonic matter is 

S dark (|z| < 1.1 kpc) < 46 M /pc 2 . (3) 

For the distribution of matter within the disk we use the isothermal sheet model (see e.g. section 
11.1 of Ref. |T 

volume distribution of DDDM as 



If the total mass of DDDM in the galactic disk is eM^y, we approximate the 



eM gal 

p(R, z) = — ex V (-R/R d )sech 2 (z/2z d ). (4) 
8nR d z d 

Here z parameterizes height above the midplane of the disk, while R is the radial direction within 
the disk. We assume the DDDM disk has a scale radius comparable to that for baryons, R d ~ 3 
kpc |34| . The value of R relevant for the measurements is the distance of the Sun from the galactic 
center, about 8 kpc. We will discuss the expected values of the disk scale height z d in Section [6] 
For now, we only need to assume z d <C 1.1 kpc, in which case the surface density does not depend 
on z d : 

eM gal 

£ d isk(M < 1.1 kpc) = — exp{-R/R d ). (5) 

lnR d 

Given this functional form, from the the surface density bound, Eq. [3] is a constraint on the fraction 
of all the dark matter that is allowed to be in a thin disk: 

e < 0.05. (6) 

This is a key result of our paper: the mass of the DDDM disk can be on the order of five percent of 
the total mass of the Milky Way. Up to order-one uncertainties, this means the mass of the DDDM 
disk can be as large as the mass of the baryonic disk, and that DDDM can carry comparable 
energy density to ordinary baryonic matter. It will be very interesting to explore whether improved 
measurements could detect new structures like DDDM disks. For instance, the ambitious plans of 
the Gaia satellite (see |35| and references therein) to produce an extensive map of the kinematics 
of a billion objects in the Milky Way could lead to a powerful probe of dark structures within the 
galaxy. 
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Other bounds can in principle arise from bounds on compact objects. Once a sufficiently cold 
disk has formed, further structure can develop within the disk. Depending on details of DDDM 
chemistry and molecular cooling that are difficult to calculate, these structures could range from 
large gaseous clouds down to "DDDM non-nuclear-burning stars" that radiate dark photons as the 
matter within them annihilates. We can estimate the size of large clouds that form within the cold 
disk based on the Jeans mass, where we treat the clouds as a monatomic gas with sound speed 



vr / 5ttT \ 3/2 /1\ 1/2 



Of course, once clouds above the Jeans mass begin to collapse, atomic and molecular cooling pro- 
cesses could lead to formation of much smaller structures. 

We summarize existing bounds on MACHOs (Massive Compact Halo Objects). For the largest 
structures above 1O 6 M0, constraints arise from heating of the disk by gravitational scattering of 



stars on MACHOs 36 . For structures above about 100 M@, including Jeans-scale clouds, the best 
constraints arise because MACHOs could disrupt wide binary star systems [37). Smaller objects 
below 100 Mq are constrained by microlensing surveys such as MACHO, EROS, and OGLE. These 
surveys are reviewed in Ref. [38] and some of their implications for dark matter are discussed in 
Ref. [39]. Observations have looked toward the Magellanic clouds, which are at relatively high 
galactic latitude and not likely to constrain DDDM. Other observations toward the galactic center 
could be more interesting. In the case of ordinary dark matter distributed throughout the halo, for 
a wide range of masses the MACHO bounds exclude the possibility that more than 10 to 15% of 
the halo consists of objects of a given mass [37]. However, we have to be careful when applying 
these bounds to DDDM, since compact DDDM objects are localized in a disk that can be thinner 
than the baryonic disk and furthermore is unlikely to be precisely aligned. Therefore the limits on 
larger objects that derive from the interaction of baryonic matter and dark matter, even if purely 
gravitational interactions are assumed, need not apply. For objects smaller than the DDDM Jeans 
mass, we clearly cannot say anything without a more detailed understanding of substructure. 

We conclude that while MACHO-type bounds might ultimately detect or constrain our scenario, 
there is no hard limit at present. The Oort limit is easier to interpret as a bound on DDDM, and 
we take it as our sole constraint on e. Perhaps future work on compact DDDM objects, together 
with new analyses of existing data, could lead to stronger bounds or detection prospects. 



3 Early Thermal History (Before Structure Formation) 

For the purposes of this paper, we consider a simple PIDM model, with a new abelian gauge group 
U(1).d, with fine structure constant ocd, interacting with two matter fields: a heavy fermion X and 
a light fermion C (for "coolant," as we will see in Section [5J, with opposite charges qx = +1 and 
qc = —1 under U(1)d. Of course, one could also consider the case that X and C are scalars. An 
interesting generalization that we will discuss briefly is the possibility that the new gauge group is 
nonabelian, SXJ (N)d, with X in the fundamental and C in the antifundamental representation. In 
the nonabelian case we make a self-consistent assumption that the confinement scale is far below the 



temperatures relevant for the phenomena we study. (In the SU(2) case, one can introduce a global 
symmetry to play an analogous role to the distinction between fundamental and antifundamental 
representations . ) 

In this section we will discuss the thermal history of the dark sector, including the amount of 
dark radiation and the abundance of X, X and C, C particles. The result is that the predicted dark 
radiation is allowed by current bounds on the number of relativistic degrees of freedom at the time 
of BBN and of the CMB, but large enough that Planck can see an interesting signal. 

We will also show that the thermal relic abundance of X and X can be of the order of the Oort 
limit, comparable to baryon density. This relic symmetric population of X and X can annihilate 
and provide an indirect detection signal. However, the light particles C and C annihilate away 
efficiently in the early universe. We therefore have the additional requirement that there is a non- 
thermal asymmetric abundance of X and C that survives to late times, analogous to the nonthermal 
abundance of protons and electrons in the SM. 



3.1 The Temperature of the Dark Sector 

The light degrees of freedom in our scenario introduce constraints (or possible signatures) since at 
early times they were relativistic and affected the expansion rate of the universe. At the time of 
BBN, the thermal bath of dark photons and also of the light species C, C will add to the total 
amount of relativistic energy density. At the time of last scattering in the visible sector, only 
the dark photons will be relativistic. The bounds from BBN and the CMB on relativistic degrees 
of freedom are usually phrased in terms of the number of effective neutrino species, so we will 
now calculate the expected number of effective neutrino species present in our model, assuming a 
sufficiently high decoupling temperature that we justify in Appendix |A| 

Suppose that, at early times, the DDDM sector and the Standard Model were in thermal equi- 
librium. After decoupling, the entropy density should be separately conserved in the visible and 
dark sectors. This means that 

„dec dec 

"*s,D _ y*s,vis /q \ 

g^MMW ~ s«,vis(t) 

with £ = (Td/T v ; s ) and g* a the effective number of degrees of freedom contributing to entropy 
density. The subscript D refers to dark sector degrees of freedom. Note that £ is, in general, a 
time-dependent quantity, as (for example) the visible sector temperature will increase relative to 
the dark sector temperature whenever visible degrees of freedom decouple from the thermal bath. 
Suppose that decoupling of the hidden and visible sectors occurs at temperatures below the W mass 
but above the 6-quark mass, which is the case if all the mediator particles have weak-scale masses. 
At this time, c vis = 86.25. The dark plasma, at this time, will contain the dark photons and C, C 
particles, leading to g^ c D = 2+ | x 4 = 5.5. It is also interesting to consider the generalization to an 
SU(iV) dark sector with C in the fundamental representation, for which g^. c D {N) = 2 (N 2 — l) + XjV. 
In the visible sector, at the time of BBN we take gf^jj = 10.75, while we expect the dark sector 
degrees of freedom to be unchanged. This leads to 

/in 75\ x / 3 
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The number of additional effective neutrino species is determined by (^d^^bbn) = 
leading to: 

AJV^n = 0.20 for U(1) D and 

AN^ = 0.077V 2 + 0.12iV- 0.07 for SU(N) D . (10) 

Numerically, AN^ is 0.46 in the SU(2) D model, 0.94 in the SU(3) D model, and 1.56 in the 
SU(4)£> model. Ref. |40| derives a conservative bound on extra-degrees of freedom during BBN, 

AN e™ < 1-44 at 95% C.L., (11) 

so the U(l)x) model is easily safe. The S\J(N)d model satisfies the bound for N < 4, with N = 4 
barely outside the 95% confidence region but easily inside if we assume decoupling at temperatures 
above the top quark mass when g^g C vis = 106.75. For an alternative point of view, we can relax 
our assumption about the decoupling temperature and ask: for what value of c vis is the BBN 
constraint saturated? It turns out that as long as 

S^vis > 19-3 (12) 

the bound is satisfied for the abelian model. This is the number of degrees of freedom when 
T^f » 200 MeV. 

An equally significant bound on the number of radiation degrees of freedom comes from the 
CMB. A recent analysis of nine years of WMAP data |4l] combined with the terrestrial experiments 
SPT |42| and ACT |43| and baryonic acoustic oscillations constrains 

an cmb < L6 at 95% Q l (13) 

At the time of last scattering in the visible sector, we have g^^^ s = 3.36 (from photons and the 
colder neutrinos) and g*^p = 2 (from dark photons) or 2(N 2 — 1) (in the nonabelian case). At this 
time the temperature ratio is 

« - (¥ ^r^™'- 

« = ( 2(N ]-?\i N X **) VS for Sm)D , (14) 
s I 2(iV 2 -l) 86.25^/ v ' y ' 

Robustly, if the two sectors are in thermal equilibrium near the weak scale, we expect the dark 
photon temperature to be around half the visible photon temperature. Alternative cosmologies, for 
instance with decoupling at much higher temperatures below which many new visible-sector degrees 
of freedom exist, could allow much smaller £, but we will generally take £ ~ 0.5 throughout the 
paper. 

The temperature of dark recombination (formation of dark atoms from X and C ions) is about 
a factor of ten below the binding energy Bxc ~ o? D vtic- Large a/j suppresses the thermal relic 
abundance of X, X and larger mc prevents efficient cooling, as we will see in Section [H] Hence, we 
favor parameter space at small Bxc where recombination in the dark sector doesn't happen until 
after last scattering in the visible sector. This means that when the CMB is formed, dark photons 
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are interacting with the dark fluid of X and C particles, with a speed of sound slightly less than, 
but of order, c/y/3. Much as for ordinary baryons, there will be dark acoustic oscillations and other 
effects from this nontrivial coupling of radiation to matter. Although not entirely correct since the 
additional degrees of freedom in our model are not yet free streaming, we interpret the bound on 
the number of effective neutrino species as a bound on free dark photons, ignoring the coupling to 
the fluid. We expect that, because the sound speed is of order the speed of light, this will be a 
good approximate guideline to whether the theory is allowed by the current data. It will be very 
interesting to do a more careful analysis that can distinguish this scenario. 

The number of additional effective neutrino species is determined by g* s .D^{tcMB) — (n)^ 3 x 
| x 2 x AN™*, leading to: 

AiV c c ff MB = 0.22 for U(1) D , 

AN^ = 4.4(iV 2 - 1)£ 4 for SU(iV) D . (15) 

Numerically, AA^ B is 0.49 in the SU(2) D model, 0.91 in the SU(3) D model, and 1.45 in the 
SU(4)£> model, so the bound is satisfied for N < 5. In the abelian model, if we ask what value of 
9*s vis saturates the CMB bound, we find it is 

S* d s e , C vis > 19.7, (16) 

essentially the same bound we derived from BBN. Thus, the abelian model is allowed provided the 
two sectors decoupled at temperatures above the QCD phase transition. The S\J(N)d is allowed 
for N < 4, and predicts sizable deviations in the number of effective neutrino species. The Planck 
satellite is expected to measure AN^^ with a one-sigma error of ±0.26 44-46 . This means 
that very soon, Planck can rule in or out SXJ(N)d models decoupling from the SM at the weak 
scale, especially for larger values of N. Improved measurements of AN^^ will also come from 
ACTpol (47] and SPTpol (48). 

3.2 Relic Abundance of X and C 

Having considered the relic radiation, we now consider the relic abundances of X and C. The thermal 



relic abundance of particles charged under a hidden U(1)d has been discussed in Refs. (13-15 
Depending on whether the mediator particles coupling X to the Standard Model thermal bath are 
heavier or lighter than X, the dark sector may be at precisely the same temperature as the SM 
when X freezes out, or as we saw in the previous subsection it could have about half the Standard 
Model temperature if the two sectors have decoupled. In Figure [T] we have plotted the curve in the 
(mx, ckd) plane which predicts e = 0.05 for the thermal relic abundance of X and X, assuming the 
SM and the hidden sector are still at the same temperature at the time of X decoupling. Taking 
the hidden sector to be at half the SM temperature leads to slightly lower values of ap. Values of 
old below the line in Fig. [I] lead to a large relic abundance that violates the Oort limit discussed in 
Section [2] whereas larger values of old are allowed only with a nonthermal mechanism for generating 
more X particles. 

Figure [l] clearly shows that we can achieve a relic density of X and X particles that saturates the 
Oort bound for reasonable values of the coupling old- It is possible that the relic density estimated at 
freezeout is later decreased by two processes with similar rates: Sommerfeld-enhanced annihilation 
of X and X at low temperatures and recombination into XX bound states which annihilate away. 
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Figure 1: ajj that yields a thermal relic abundance of X, X that is a 5% fraction of the total DM density 
for different mx- 



Using the recombination rate given in 15 49 , we plot in Fig. [2] the curves Tree — H for particular 
choices of visible and dark sector temperatures. These show that recombination and subsequent 
annihilation of XX bound states does not wash out the abundance of X, X for < 0.01 and 
mx ~ 1 GeV. 

The light species C with mc <C mx freezes out at much later times, and has a much larger 
annihilation rate than the heavy species, by a factor (mx /mc) 2 ■ As a result, the thermal relic 
number density of C is much smaller than that of X, by a factor mc/mx- This means that we 
expect any symmetric component of C and C to annihilate away almost completely at dark sector 
temperatures a factor of 20 below the C mass. The existence of light C particles is crucial to 
dissipative dynamics, as we will see in detail in Section [5] This means that only a nonthermal 
mechanism for producing C particles can be consistent with dissipative dynamics. 

We assume asymmetric nonthermal abundances of DDDM, in which we have a net C number 
n C ~ n c 7^ 0- The universe should be charge symmetric and this means that there should also exist 
an asymmetry in X, nx — n x = nc — n c - We assume that n x = n c = at late times, whereas X 
and C survive. The idea of asymmetric dark matter has inspired many proposals for how the dark 
matter and baryon abundances may be related 1 50- 60 1 , and similar considerations could explain the 
amount of DDDM. We expect for reasonable parameters that the asymmetric abundance of X and 
C will be supplemented with a relic symmetric population of X and X, allowing for the prospect 
of interesting indirect detection signals when X and X annihilate to Standard Model particles. 

It is possible that other nonthermal scenarios could leave a symmetric population of C and C that 
survives to late times. For instance, late-decaying dark sector particles could produce additional 
C particles after thermal freezeout of C. However, even these late-produced C particles could 
annihilate at temperatures below B cc = a 2 D mc/^, through Sommerfeld-enhanced annihilation or 
recombination. One way to prevent this would be if they were first bound up into dark atoms 
with X, which have a larger binding energy and might protect the C particles from annihilation. 
However, recombination of XC dark atoms is much slower than naively expected, because one 
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Figure 2: Above the curves, the recombination rates are larger than the Hubble rate, leading to XX 
annihilation that depletes the abundance of the symmetric relic component. 



dark atom can be ionized by the dark photon emitted when another dark atom is formed. This is 
analogous to what happened for ordinary hydrogen in our universe [61]. It has been studied for 



dark atoms in Refs. 16, 18 . The result is that CC annihilation would always be much faster than 
the sequestering of C particles inside dark atoms with X. 

A loophole arises in the case of nonabelian gauge theories. In that case, the dark gluon emitted 
by formation of one dark atom could ionize another dark atom, but it could also first encounter 
another dark gluon and scatter off it. In the process, one of the two dark gluons loses energy. 
Because the number of formed atoms and the number of dark gluons emitted in recombination 
must be equal, we expect an order-one fraction of dark gluons can lose energy and thus become too 
soft to ionize a dark atom. This could allow the C and C particles to be stored in tightly-bound XC 
and XC bound states, possibly allowing an interesting symmetric abundance of C, C to survive to 
late times and contribute to cooling. A full numerical analysis of this rather complicated cosmology 
is beyond the scope of this paper. 

Because the cosmology in which we simply assume an asymmetric component of X and C is 
simpler, through most of the paper we will assume asymmetric dark matter. 



4 The Fermi Line and Other Indirect Detection Signals 

The DDDM scenario was motivated in part by the observation of gamma-ray line emission at about 
130 GeV (135 GeV after energy recalibration) in the galactic center using 3.7 years of Fermi-LAT 
data |3|El|62] . The observation is not decisive, and certain features in the data suggest it could very 



well be an instrumental effect [63-66 . Nonetheless, the suggestion serves as a concrete example 
where a possible enhancement of an indirect signal inherent in our scenario could be critical. In 
fact, any high-energy gamma ray line observable in the near future and consistent with continuum 
bounds would require similar phenomenology. Furthermore, a large flux of high-energy electrons 



and positrons recently observed in cosmic rays by PAMELA and subsequent experiments 67-71 



also requires large boost factors to explain in terms of dark matter annihilation [72] (though it is 
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plausibly due to pulsars). This shows that, even apart from the specific case of the Fermi-LAT line, 
it is very interesting to consider general mechanisms that can produce large boost factors. We will 
see in this paper that boost factors as large as 10,000 are conceivable for DM of mass 130 GeV 
with U(l)x) coupling in the range to give the right thermal relic density. The large astrophysical 
enhancement we find will be due to a very thin dark matter disk. In modified models, related 
enhancements could explain boost factors needed for PAMELA or a possible AMS signal |73| as 
well. 

We first consider the required signal enhancement in the context of a simple model in which 
DM, (f) and the new charged particles S are all scalars. 

-C D X^lsf + mllSf + XslSf + mll^ + X^, 

in which <p is charged only under U(l)£> and S has charge 1 under the usual U(1)em- One could 
also consider models with fermionic DM. In that case, to avoid kinetic mixing between U(1)d and 



U(1)em 74 , one could design a resonant annihilation model where DM annihilates through an 
intermediate boson and a U(1)em charged particle loop to two photons, although anomaly-like 
constraints on charge assignments can make such models safe (as explained in Appendix [A]) even 
when a particle is charged under both U(l)s. 

The observed photon line could be consistent with DM particles annihilating at one loop to 77 
with an unexpectedly large cross section of order (av) ~ 10~ 27 cm 3 /s. In our example model, the 
cross section of the DM annihilation to diphotons would be 



B 



™^-m 327r 3 m 2 



aA</,5T A)(T</>) , (17) 



where B is the boost factor which could either come from microscopic physics such as Sommer- 
feld enhancement, which we will discuss in this section, or from astrophysics, such as the density 
enhancement we will discuss in the following sections. Here 

A (t) = -t + t 2 /^- 1 ) with T 4> = m\ /m% (18) 
f(x) = arcsin 2 -v/x. (19) 



(For ms < m^, it is necessary to analytically continue f(x).) Demanding (av)((f>^(f> — > 77) = 
10~ 27 cm 3 s _1 and fixing A^s = 1, one can derive the required B for a given ms- The result is 
presented in Fig. [3] Thus to explain the line, without relying on large couplings and tuning ms to 
be close to the DM mass, one needs a huge boost factor. It is also difficult to make this scenario 
consistent with a thermal relic abundance. The observed av to photons is too small, leading to 
an overabundance of <f) if this is the only annihilation channel. But adding larger annihilation 
channels, such as — > W + W~ , is typically in tension with the absence of an observed gamma 
ray continuum. One possible explanation is if S is slightly heavier than (/), so that (fxfy — > SS^ can be 
an important annihilation mode in the thermal environment of the early universe, but is impossible 
today |6 75 



For DM charged under an unbroken \J(l)n, there is indeed an enhancement, the Sommerfeld 
enhancement, 

_ 2ira D Q 2 D /v 
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Figure 3: Fixing (av)(4>4> — > 77) = 1CP 27 em's' 1 . X^s = 1 and = 130 GeV, the boost factor B needed 
to generate a signal that current Fermi line search is sensitive to. 



where v is the DM velocity. However for a D Q 2 D < 0.1 and v 1(T 3 , S < 1000. If one fixes the DM 
thermal relic J7/i 2 to be 0.11, ar>Q 2 D ~ 3 x 10~ 3 , S ~ 20. Thus Sommerfeld enhancement itself is 
not sufficient enough to get the desired annihilation cross section. 

In what follows, we will see that enhanced density from DDDM could be sufficient to generate 
such a large boost factor. Of course the precise value for the boost will depend on the precise 
parameters of the dark matter candidate (and some as-yet-unknown astrophysics) as we discuss 
below. 



5 Cooling 

The enhanced signals we discuss arise as the result of the interacting component of dark matter 
collapsing into a disk. We now consider when and how this can occur. This is equivalent to the 
question of when cooling is sufficiently quick to allow for collapse, so we now investigate the question 
of how interacting dark matter can cool. 

The cooling has many features in common with ordinary baryonic matter. DDDM first adi- 
abatically cools through the expansion of the universe. As with baryons during the formation of 
a galaxy, the interacting dark matter will already be present (in the primordial overdense region 
that seeds an early galaxy halo or in progenitor halos that merge into a larger galaxy) and will also 
accrete onto the galaxy from the intergalactic medium. After virialization through shock heating, 
baryons cool from different processes: atomic and molecular interactions, Compton cooling, and 
bremsstrahlung radiation. All of these require light electrons in order to have a sufficiently rapid 
rate. 

The same mechanisms will be required for DDDM: cooling occurs sufficiently rapidly only when 
a light particle is present that also interacts under the dark U(l) (or more generally, whatever force 
is relevant). Therefore at the time of the initial accretion, part of the DDDM might be bound into 
atomic-like states of heavy and light dark matter. As discussed in the previous section, for instance, 
the simplest model is an asymmetric population of X and C, which like ordinary hydrogen can form 
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bound states in the early universe with some residual ionization (as calculated carefully in [18]). A 
relic population of X and X may also survive, so the initial conditions will involve a mix of dark 
atoms and dark ions. However, as we will now see that shock heating will destroy any initially 
bound atoms, we can consider cooling in this section without determining the exact fraction of 
bound states in the very early universe. 

As dark atoms fall into the overdense region, their gravitational potential energy converts to 
kinetic energy. Initially they are quite cold, but when falling into the galactic center, particles slow 
down as they encounter other infalling particles, forming a shock wave which expands outward, 



containing pressure-supported gas inside 76 77 . This shock-heating process converts the kinetic 



energy of the DDDM gas to thermal energy at the virial temperature, 

G N Mfi o £ i v M V 110 k P c <n\ 
= W W 8 ' 6 keV ^ ToOGeV^— • (21) 

where M stands for the mass of the virial cluster and /i = p/n\s the average mass of a particle in the 
DDDM gas. We have taken a fiducial value for the mass of dark matter in the Milky Way galaxy, 
= 1O 12 M . This is reasonable since the initial density perturbation induces gravitational 
collapse in the dominant dark matter component for which neither baryons nor the subdominant 
interacting dark matter should be very relevant. Note that for a virial cluster of the same mass 
and radius, DDDM will be much hotter than baryonic matter, with a temperature enhanced by 
~ mx/rrip- The binding energy of the ground state of the dark atom is 

Bxc - ^p, (22) 

less than or of order the binding energy of ordinary hydrogen, so we expect T v ; r 3> Bxc- At these 
temperatures the DDDM in the virial cluster will be completely ionized, even if it had recombined 
into dark atoms or dark molecules before virialization. Hence we can start off thinking of free X 
and C particles. 

The same cooling processes that apply to baryons potentially apply to DDDM. An ionized 
dark plasma in the virial cluster can be cooled through bremsstrahlung and Compton scattering 
off background dark photons. Compton scattering is more efficient at larger redshift, when the 



dark photon background was hotter. Based on the results of Section 3.1, we take the dark photon 
temperature to be Tp ~ 0.57cmb- (This is the temperature of the dark cosmic background photons, 
which is to be distinguished from the temperature of X and C particles in the galaxy, T va .) The 
timescale of the bremsstrahlung cooling is 



o , 3/2^1/2 

^brem ~ ., „ q 

16 nxnc or D 



(23) 



10 4 yr 



f r vir cm 3 / «em\ 3 / mc N 
K n c \ a D J V m e , 



where in the second line, we assume iix = nc for simplicity. At the end of the section, we will relax 
this assumption. This time should be compared to the age of the universe in order to show that the 



14 



cluster efficiently cools down. The timescale for cooling through Compton scattering is 



_ 135 n x + n c m% 
C ° mpt0n ~ 64.3 nc a 2 (T o (1 + 2)) 4 W 

« 4xl0 1 V n * +reC f^W „ 2K r Vf^V, 

n c V «o y v^(i + z)y V"W 

where T D is the current dark CMB temperature and z is the redshift. In Figure gat left, we show 
contours in the plane of mc and redshift along which the bremsstrahlung and Compton cooling 
rates are equal, for different choices of old- Because the Milky Way galaxy was starting to form 
before z = 2, Compton cooling of DDDM would be important within the Milky Way at early times. 
Compton scattering could also be important for smaller old and mc- We illustrate this in the 
right-hand plot of Fig. [4j which shows the contours in the (mc,ap) plane along which the two 
rates are equal and along which the faster rate equals the age of the universe. As the dark photon 
background cooled, bremsstrahlung would have become increasingly important. We use the generic 
term i coo i for whichever time scale is shorter: i coo i = min(tbremj ^Compton)- 




Figure 4: Comparison of the rates of bremsstrahlung and Compton cooling. At left: the value of mc for 
which the rates are equal, as a function of redshift. To the right of the curves, i.e. at early times, Compton 
cooling dominates. At right: the contour in the (mc, cld) plane along which the bremsstrahlung cooling rate 
equals the Compton cooling rate (black dashed line) and the contour along which the cooling rate equals 
the age of the universe (solid purple line) . This shows that Compton cooling is the dominant effect at small 
mc and ajj, while bremsstrahlung dominates for larger values. In both plots, we have taken an NFW virial 
cluster of radius 20 kpc. 

In order to verify that bremsstrahlung or Compton scattering leads to cooling, we first make 
some consistency checks. The emitted dark photons must escape from the galaxy and carry away 
energy without being reabsorbed. The primary process by which a dark photon would interact is 
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through scattering with a light C particle, so we can approximate the photon's mean free path by 



3lJlc ~ 1.5 x 10 s kpc, (25) 



ottlc 8Tia D nc 

where we have used the Thomson cross section for jd~C scattering with au = a, rac = m e , and 
mx = 100 GeV while assuming equal X and C number densities at e = 0.05 and a virial radius of 
110 kpc, namely 

n x = n c « 3.3 x lO" 6 cm" 3 ( W ^ Ge l) . (26 ) 

V m x J 

The long mean free path shows that photons readily escape the galaxy at early times. Furthermore, 
because i ~ 10 6 i? v i r , photons will continue to escape even if the initial DDDM distribution collapses 
by a factor of 10 18 in volume. This is sufficient to allow a disk to form, especially considering that 
once the DDDM assumes a disk-like shape, photons can escape more efficiently through the thin 
direction of the disk. 

We also need to check that both light and heavy particles could cool. When light particles 
scatter on heavy particles and emit bremsstrahlung photons, it is mostly the light particles that 
lose energy. Similarly, Compton scattering is dominantly scattering of the light particles on dark 
background photons. However, if heavy and light particles remain thermally coupled, the cooling 
of the light particles is sufficient. 

Thermal coupling occurs when the rate for Rutherford scattering of the light particles on the 
heavy particles exceeds the cooling rate. In this case, the heavy particles cool adiabatically, with 



scattering keeping the light and heavy species in in kinetic equilibrium [78 , 79 . The timescale for 
this equilibration process is 



m x m c (E c I mcf 12 (0 „, 



D n c log 1 + 



2/3 

2 , .5 



4.3xl0 4 yr — _^L- 2 ^1 

\a D J VlGeW \m C ) n c 
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ct \ / mx \ 2 ( rn e \ 2 cm 10 



log 1 + ^ 

a D n C 



where Eq is the kinetic energy of the light species; in the second line, we take Ec/mc = 3T v i r /mc- 
In part of our parameter space, t eq <C t C ool and the light and heavy species cool adiabatically 
together. 

Rutherford scattering has a enhancement, but when mx/mc is very large, in a thermal 

system vc is not small. Thus for large mx/mc, as well as in the region of parameter space where 
old is very small, the equipartition time from two-body scattering processes is not sufficient to cool 
the heavy particles. In this case, we expect cooling should still occur but that cooling involves 
nonequilibrium physics, at least initially. If the light particles contract as they cool, while the heavy 
particles are unaffected, a charge separation would occur between the larger cloud of X particles and 
a smaller cloud of C particles. This wwould produce dark electric fields that pull the X particles in. 
It would be interesting to simulate or model more completely the resulting dynamics, but it seems 
inevitable that, since cooling continues to rob the system of kinetic energy, eventually both X and 
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C will cool. As they contract into smaller volume, larger values of nx,c make Rutherford scattering 
more efficient, and the cooling process will eventually be describable again by equilibrium physics. 

Hence, we work under the hypothesis that whenever the cooling time scale i coo i is less than 
the age of the universe, cooling occurs. At this point we should mention one further subtlety: 
equipartition will speed up the light particles relative to the heavy ones by a factor \fm~x~jm~c, 
and so for sufficiently small mp we should use the formula for relativistic bremsstrahlung rather 



than Eq. 24 Since the rate of energy loss from relativistic bremsstrahlung exceeds that from 



nonrelativistic bremsstrahlung by a Lorentz factor, cooling will become only faster |49|. Thus 



Eq. 24 is a conservative estimate. 

In Fig [5] we present contours on which the cooling timescale is sufficiently rapid. We derive the 
bounds by assuming two different number densities tiq = Tlx- First we assume conservatively that 



DDDM is uniformly distributed over a 110 kpc sphere as in Eq. 26 The bounds for mx = 100, 1 
GeV are shown in the upper row of Fig. [5] Because DDDM will tend to fall into the halo, we expect 
that in fact cooling will be more rapid due to enhanced number density in the halo's central region. 
To obtain a more optimistic estimate, then, we estimate the time scale again using an NFW profile 
with a characteristic scale R s = 20 kpc. We still use the virial theorem 

13Q„M(fl.) = 3 

2 5R S 2 ° V ' 

where T20 is the temperature inside this region with radius R s = 20 kpc and M(R S ) is the mass 
inside the region. Then the number density is chosen to be 

3M(R S ) 1 . , 

" c = "- = e W^- (29) 

The results are shown in the bottom row of Fig. [5] and potentially allow masses an order of magni- 
tude larger than the conservative estimate when bremsstrahlung is the dominant cooling mechanism. 
These plots show when efficient bremsstrahlung and/or Compton cooling can begin. Once cooling 
begins, increased density makes it more efficient, so the process will continue. For smaller mc and 
old, Compton cooling could be faster than bremsstrahlung cooling. Thus the curves in Fig. [5] have 
kinks which correspond to transitions from bremsstrahlung cooling domination to Compton cooling 
domination. We expect that both bremsstrahlung and Compton cooling will continue until heavy 
and light ions become cold enough to recombine into dark atoms. 

From the left-hand plots in Fig. [5] we see that for mx = 100 GeV, there is a small region 
of parameter space where mc ~ 1 MeV and old ~ 0.1 where bremsstrahlung cooling happens 
within the age of the universe and the X and C particles are in equilibrium. A much larger region 
of parameter space with smaller mc has a fast cooling time but slow equipartition, so a better 
understanding of nonequilibrium cooling is needed to be certain of the fate of DDDM in this region. 
On the other hand, the right-hand plots of Fig. [5] show that a much larger region of parameter space 
cools adiabatically when mx = 1 GeV. In particular, a Standard Model-like choice mx = 1 GeV, 
mc ~ rn e , and old ~ ol is on the edge of the region that cools efficiently by bremsstrahlung. (The 
SM gets some help from the baryon abundance being 15% of all matter rather than e = 0.05 in our 
plot.) A larger region of parameter space down to mc = 10 keV with old between 10 -4 and 10 _1 
can cool efficiently, typically through bremsstrahlung but with Compton cooling predominating at 
the smaller values of old- 

For the evaluations above, we assumed nc = nx, as is the case for fully asymmetric dark matter. 
This is not necessary. In Fig. § we fix otD such that the thermal relic abundance of X is 5% of the 
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Figure 5: Cooling in the {mc,OiD) plane. The purple shaded region is the allowed region that cools 
adiabatically within the age of the universe. The light blue region cools, but with heavy and light particles 
out of equilibrium. We take redshift z = 2 and Tjj = Tqmb/2- The two plots on the left are for mx = 100 
GeV; on the right, mx = 1 GeV. The upper plots are for a 110 kpc radius virial cluster; the lower plots, 
a 20 kpc NFW virial cluster. The solid purple curves show where the cooling time equals the age of the 
universe; they have a kink where Compton-dominated cooling (lower left) transitions to bremsstrahlung- 
dominated cooling (upper right). The dashed blue curve delineates fast equipartition of heavy and light 
particles. Below the dashed black curve, small an leads to a thermal relic X, X density in excess of the Oort 
limit. To the upper right of the dashed green curve, Bxc IS high enough that dark atoms are not ionized 
and bremsstrahlung and Compton cooling do not apply (but atomic processes might lead to cooling). 



18 




lO -7 10 -6 10 -5 10 -4 10 -3 10" 2 10"' 1(T 6 10 _s 10 -4 10~ 3 10 -2 



m c [GeV] m c [GeV] 

Figure 6: Fixing old to yield the thermal relic abundance of DDDM as 5% of the total DM relic abundance, 
the minimal nc/nx with a cooling time scale equal to the age of the universe as a function of mc for 
mx — 100,1 GeV. We choose redshift to be z = 2 and dark temperature To = Tcmb/2. The bounds 
are from the Compton cooling process, which for the chosen mx and ajj dominates (so the bounds are 
independent of the DDDM density profile). 



total DM density (as in Figure [T]) and plot the minimal nc / nx needed to have a cooling time scale 
shorter than the age of the universe. Combining both plots, we see that in most of the parameter 
space, sufficient cooling requires light DDDM with density greater than thermal (but comparable 



to the thermal abundance for X), which we will assume to be present. As discussed in Section 3.2 
we take the nonthermal component to be asymmetric. 

Note also that in principle other cooling processes might occur. We expect atomic or even 
molecular processes or collisional cooling won't be important until below the expected recombination 
temperature, which is below the binding energy. However, note that it is only important at this 
point to establish that cooling can indeed occur, and that the temperature can be sufficiently low 
to form a thin disk no bigger than that for the baryons. Bremsstrahlung or Compton scattering 
with sufficiently light C particles ensures this can indeed be the case. 



6 Disk Formation 

Having established that DDDM can efficiently cool via bremsstrahlung or Compton scattering, we 
now consider how it will be distributed within the Milky Way. Like any matter falling into a halo, 
DDDM will have angular momentum, and so, as with baryonic matter, we expect that it will cool 
into a rotationally-supported disk. Because DDDM does not have supernova feedback and other 
processes that may be important in the evolution of the baryonic disk, this is not entirely obvious, 
and we rely on recent sophisticated numerical simulations in which disk formation occurs without 



including stellar and supernova feedback 80,81 , rather than the earliest simulations in which 



baryons formed small clumps rather than disks [82]. It is important to have further numerical work 
to confirm that this is true; DDDM that could form clumps instead of disks could also be extremely 
interesting. 
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We assume the disk mass distribution from Eq. [4] Assuming the disk scale length Rj is much 
larger than the scale height Zd, we can neglect radial derivatives in the Jeans equation for an 
axisymmetric system and estimate, at the galactic center: 



1.2 



G N eM^ 



10 



-3 



Rd 

e 



(30) 



where v v z ms = is the velocity dispersion of DDDM in the vertical direction, p cen ter = p(0, 0) is 
the central mass density, and we have used Eq. [4] in the second step. This estimate assumes the 
gravitational potential is dominated by the disk, i.e. it ignores the effects of baryons and of ordinary 
dark matter, but this is self-consistent to the extent that the disk is quite thin and thus the DDDM 
density is locally much larger than that of baryons and ordinary dark matter. We estimate that the 
vertical velocity dispersion corresponds to the temperature at which cooling stops, v\ ~ ^cooled/ m x- 
We also assume that the disk scale length Rj is comparable to that for baryons, around 3 kpc |34| . 
Eq. 30 should be viewed as a rough estimate; in particular, the detailed spatial distribution of 
DDDM may not precisely correspond to Eq. |4j 
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Figure 7: Estimates of cooling temperature T C ooied/ Bxc in the {mc,OLd) plane. Black solid curves: ioniza- 
tion fraction x = 0.1; green dashed curves: x= 0.01. 



In order to obtain concrete numbers from Eq. 30 , we need an estimate of the final temperature 



^cooled- Bremsstrahlung and Compton cooling will cease to be efficient once the light particles 
are slow enough to recombine into dark XC atoms, at temperatures low compared to the binding 
energy -BxcQ A rough estimate of the temperature at which this occurs is found by solving the 



'We are assuming an asymmetric scenario in which X and C are present, possibly with a symmetric X and X 
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Saha equation, 



2 1 /TW\ 3 / 2 

J exp(-^j, (31) 



n x n c _ x 2 1 ( Tm c \ ( B X c \ 

nxcn 1 — x n \ 2tt 



where nxc is the bound state number density and n = nxc + n x an d nx = nc- The ionization 
fraction x = nc/n. In order to obtain the relevant density, we assume that the gas has already 
cooled into a disk with scale radius and scale height Zd 

= _Po_ _ eM gai ^ G N (eM gal ) 2 
n m x 8irR 2 d z d m x 128vri^T ' 1 ' 

Combining Eq. [31] and Eq. [32] and requiring the ionization fraction to be smaller than 1, e.g., 0.1 
or 0.01, we find the results shown in Figure [7] which we summarize as: 

T C ooied ~ (0.02 - Q.2)B XC - (33) 

Thus, we expect that cooling stops at a temperature of about 10% of the binding energy. This leads 
to estimates of the disk scale height that are substantially thinner than the baryonic disk. Over the 
bulk of parameter space, we find that the results are well-described by a power law: 

-or; f a D\ 2 m c 100 GeV 
^ 2 - 5pC U2) 10-3 GeV m x (34) 

The 1 /mx scaling arises because, at a given temperature, the velocity of the dark atoms is smaller at 
larger X masses. In other words, we expect large boost factors for weak-scale dark matter because 
it is much heavier than baryons. For lighter X particles, smaller values of atp can still allow cooling 
and a thermal abundance of the symmetric component, in which case we again can get large boost 
factors simply because the binding energy, and hence the temperature of dark atoms relative to 
baryonic atoms, can be much smaller. 

Further cooling could occur, as in the baryon sector, through molecular processes; on the other 
hand, heating processes could also occur that would thicken the disk. For instance, the gravitational 
influence of interstellar clouds on the vertical distribution of stars in the Milky Way is important |83| 
[84] , A molecular cloud that accelerates stars will also accelerate gas particles, like the XC bound 
states, and this could thicken the disk. However, stars are collisionless, while XC bound states could 
cool down again, so we expect the velocity dispersion imparted by interstellar clouds to be smaller 
for DDDM than for stars. In the absence of a more thorough treatment of such possible heating 



mechanisms, we can only say that the true disk thickness is expected to lie between Eq. 34 and 
the height of the baryonic gaseous disk (on the order of 100 parsecs). It would be very interesting 
to see if simulations could provide a more robust estimate of the disk height, which is crucial for 
understanding the possible enhancement in dark matter detection signals. 

The angle between the baryon and DDDM disks also plays a key role in the observability of 
DDDM, especially for direct and indirect detection. We would expect that gravity would tend to 
align these structures in a timescale set, very approximately, by t ~ Ryf R/GM^^ ~ 10 7 yr. In 
fact, even as the galaxy first formed, the angular momentum vectors of baryons and DDDM could 
already have been approximately aligned, because filaments in the cosmic web define preferred 



population, but C is absent. As discussed in Sec. 3.2 a symmetric scenario may be possible in models with a more 



complex cosmology, in which case XC recombination and CC annihilation could also occur during cooling. 
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directions for accretion. Recent numerical simulations of the galaxy |85] have found that the stellar 
and gaseous components of the baryonic disk are typically aligned to within about 7°, and the 
angular momentum vector of dark matter in the inner halo is somewhat less aligned, with a median 
angle of 18° to the angular momentum of the gaseous disk. The fact that simulations see a much 
better alignment of the angular momentum of the baryonic disk with the angular momentum of 



dark matter in the inner part of the halo 85 86 , rather than the entire halo, is reflective of the 



gravitational alignment that we expect to happen between the two disks in our model. Because we 
expect approximate alignment, indirect detection signals from the galactic plane might be expected 
as we discuss in the following section. 



7 Indirect and Direct Detection 
7.1 Indirect detection 

A dramatic signal of DDDM can arise from annihilation of dark matter particles with their an- 
tiparticles, e.g. annihilation of residual ionized X with X into gamma rays for example. Because 
photons travel unimpeded to us, such a signal could in principle provide a map of the dark disk 
on the sky, giving striking visual confirmation that dark matter has cooled into a structure distinct 
from a typical halo. The gamma-ray intensity in a given direction is the line-of-sight integral of the 
DM number density squared along a given direction, 

d$ 7 _ 1 {vv) 17 2 
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with: 



J= f dbdl [ ^cos&f^V , (36) 
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where p & is the normal DM density at the Sun, p = 0.3 GeV cm -3 . d@ is the distance from the 
Sun to the galactic center (GC), d® ~ 8.3 kpc. The integral is over the region of interest (roi) at 
the GC. The smallest region centered around the GC that the Fermi-LAT experiment is sensitive to 
is a 0.2° x 0.2° square due to finite angular resolution, which corresponds to a 28 pc x 28 pc region 
around the GC. Thus for a disk height Zd > 28 pc, we expect that Jdddm scales as zj 2 . Fig. [8] 
shows the local density enhancement of DDDM compared to the normal DM defined as 

Jdddm 
Jdm 

as a function of the DDDM disk height Zd, where for normal DM, we used an Einasto profile 



/OEinasto(^) = Ps exp (-(2/oe) {{r/r s ) aE - 1)), 



(37) 



with r s — 20 kpc and = 0.17. p s is fixed to achieve the correct p@. The resulting boost factor 
arises not only from the compression of the disk in the vertical direction, but also because the disk 
scale length in the radial direction is somewhat smaller than the radial spread in the distribution 
of ordinary dark matter. 

Clearly an enhanced DDDM density would be distinguishable. Even if the density enhancement 
Jdddm /Jdm integrated over the region of interest is modest, the distribution of photons within this 
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Figure 8: Local density enhancement in DDDM, as a function of disk scale height Zd, in a square region 
around the GC fixing e — 0.05 that DDDM is 5% of the total DM density. Red: region within b C 
(-1°,1°),Z C (-1°,1°). Green: region within b C (-0.1°, 0.1°), I C (-0.1°, 0.1°) (current Fermi-LAT 
angular resolution). Black: region within be (-0.01°, 0.01°), I C (-0.01°, 0.01°). 
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Figure 9: Sky maps of the photon flux shape in arbitrary units for different DM profiles. Upper: Normal 
DM with an Einasto profile. Middle: DDDM in a disk aligned with our disk. Lower: DDDM in a disk 
misaligned with our disk by 18°. The DDDM images have disk scale height Zd — 100 pc. 



region — and especially at larger distances from the galactic center — can be radically different for 
DDDM and ordinary dark matter. Some illustrations of the photon flux over the sky are shown in 
Figure [9j 

Another feature of a possible indirect detection signal from DDDM annihilation is that a larger 
Sommerfeld enhancement can arise due to the smaller velocity dispersion of DDDM. It is usually 
assumed that the DM halo is approximately isothermal, and thus the velocity distribution is mostly 
Gaussian with a dispersion ~ lCP 3 c. However, DDDM travels in circular orbits around the GC 



23 



and its velocity dispersion could be much smaller than 10~ 3 c. Specifically, its velocity dispersion 
is determined by T coo i e d through v 2 = <iT coo i c ^/mx- For example, in the parameter space where 
bremsstrahlung or Compton cooling is efficient as shown in Sec. [5] v 2 < 1CP 9 for mx = 100 GeV. 
Thus the Sommerfeld enhancement factor, which scales as 1/v, could be increased by a factor of 10 
or more compared to non-dissipative DM charged under U(l)x> with the same charge and mass. 

In summary, for indirect detection, the DDDM scenario could easily accommodate a large boost 
factor from local density enhancement in the range (10-1000) depending on the disk height. Also 
due to a smaller velocity dispersion, DDDM could have a larger Sommerfeld enhancement, giving 
rise to another boost factor of 0(100). Thus the DDDM scenario could easily explain the suggestion 
of a Fermi photon line at around 130 GeV without large couplings or tuned masses. Again, we would 
like to emphasize that this photon line only serves as an example that the DDDM scenario could 
lead to interesting and distinctive indirect detection signals. 

Throughout this discussion we have assumed that the DDDM disk and the ordinary dark matter 



are centered on the same location. The Fermi 130 GeV line is arguably off center 62 , which has 
provoked some debate, with numerical simulations showing that dark matter may be displaced from 
the galactic center [87] and others arguing that tidal disruption prohibits such a displacement |88| . 
It would be interesting to explore the similar question of whether the DDDM disk and the baryonic 
disk can be centered on different locations. 

7.2 Direct detection 

Direct detection of dark matter could in principle be possible if the Earth is located within the 
DDDM disk. In the most optimistic case when the Sun is in the DDDM disk, the DDDM density 
at the position of the Sun could be as large as 6 GeV/cm 3 , 20 times as large as the normal DM 
density, for e = 0.05 and z,i = 100 pc. However, the spectrum of DDDM scattering off nucleons 
would be very different from that of an ordinary WIMP with the same mass, and we will see that the 
kinetic energy of DDDM is too low to produce a measurable signal in conventional direct detection 
experiments. 

The kinematics of direct detection involves a dark matter particle moving with a nonrelativistic 
velocity v in the lab frame, which scatters off a stationary nucleus. Depending on the scattering 
angle, the recoil energy imparted to the nucleus can take any value between zero and 



m N 

n , fiN \ 2 100GeV ( v x \ 2 
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where vx is the dark matter velocity, mx is the mass of the target atom, and fj,x is the reduced 
mass of the DDDM-nucleus system. Most experiments are sensitive to energies above a threshold 
value of En, below which noise and various backgrounds can overwhelm any possible dark matter 
signal. Having a threshold Er corresponds to being sensitive to a minimum value of the dark 
matter velocity vx- For ordinary dark matter, there is a broad spectrum of velocities that can be 
approximately modeled as an isothermal distribution with typical velocity 10 _3 c. But DDDM is not 
ordinary dark matter: after cooling, it is in the form of a rotationally supported disk, and a typical 
DDDM particle will move in a circular orbit around the center of the galaxy. Near the Sun, both 
DDDM and the solar system would be in approximately the same circular orbit, so the large radial 
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component of their velocity will be identical. Only deviations from this typical circular velocity can 
contribute to scattering. The rate for spin-independent elastic scattering is: 



dT = N '-r,~~n X*F{E R fE R e{v x - (39) 
dE R 2m x n-iv x 

where Nf, and A are the number, mass, and atomic number of the target atoms; mx, px and 
vx are the the mass, local density, and velocity of DDDM at the Sun; a n is the zero-momentum 
spin-independent DDDM-nucleon scattering cross section; p n is the reduced mass of the DDDM- 
nucleon system; F(E R ) 2 is the nuclear form factor; and u m ; n the minimum DDDM velocity needed 
to create a recoil with recoil energy E R . Before taking into account the nuclear form factor F(E R ) 2 , 
the spectrum would be flat between and E^ ax . However, the nuclear form factor F 2 (E R ) is in 
general an exponentially falling function, which suppresses higher energy recoils, yielding a falling 
spectrum with an end point at E^ ax . 

The typical threshold for current direct detection experiments is a few keV in nuclear recoil 



(e.g., the CoGeNT threshold is 0.5 keVee ~ 2 keVnr). Equation 38 shows that the velocities in the 
lab frame need to be larger than 10 _4 c in order for DDDM to produce hard enough recoils to be 
detected at such experiments. There are several sources of relative velocity between the detector 
and the DDDM. One is the peculiar velocity of the Sun, which does not have a perfectly circular 
orbit. Another is the motion of the Earth around the Sun. Both of these velocities are around 
10 -4 c, too small to give an easily detectable signal. Other sources of relative velocity could arise 
from inhomogeneities in the disk that lead to deviations from perfectly circular orbits. For instance, 
the spiral arms of the baryonic disk are density waves, analogous to traffic jams, across which the 
radial velocities of stars vary. The Sun is in such a spiral arm, and the DDDM disk may also have 
density waves or other structures in which velocities differ. However, given that the spiral arm in 
our disk only modifies stellar velocities by ~ 10~ 4 c, it is not clear that such effects can be large 
enough to change our conclusion. 

In general, it is interesting that due to the DDDM's small velocity, only the energy bins close 
to an experimental threshold could be sensitive to DDDM scattering. So far, the importance of 
energy calibration around the threshold has been mostly emphasized for ruling in or out the light 
DM scenario. Yet from the discussions above, pushing the thresholds of direct direction lower could 
also be important for the DDDM scenario, or in general, DM with a small velocity. 

Another possible way to detect DDDM directly is to look for single- or few-electron events if 
it scatters with electrons, causing ionization of atoms in a detector target material. In particular, 



dual-phase liquid xenon detectors could have sensitivity to such small ionization signals 89-91 . 

The low relative velocity would tend to suppress direct detection even in the most optimistic 
case, when the Earth is directly inside the DDDM disk. Another possible suppression mechanism 
is that the Earth could be located in a region of low DDDM density. The Earth sits about 10 to 20 
parsecs above the galactic plane |92| and about 8 kpc from the galactic center. A 5° inclination of a 
Zd = 100 pc DDDM disk would suppress the local DDDM density by a factor of 10 compared to the 
normal DM density at the Sun, p = 0.3 GeV/cm 3 . Hence, even improved low-threshold experiments 
could not completely rule out DDDM, since it is always possible that the density near the Earth is 
simply too small to observe. 



25 



7.3 Solar capture 

Another possibility for detection arises from solar capture. As dark matter particles pass through 
the Sun, they could scatter off nuclei inside the Sun and become gravitationally bound. With 
subsequent scattering (between themselves and nuclei), they could eventually accumulate in the 
center of the Sun. Captured X and X particles could subsequently annihilate into various SM final 
states. For instance, they could annihilate into Z7 or ZZ through the same loop of charged particles 
that leads to monochromatic photon lines. Zs would subsequently decay into energetic neutrinos, 
which could be observed in neutrino telescopes on Earth, such as IceCube. 

Currently IceCube constrains the spin-independent capture rate to be C|J ~ 10 22 s~ 4 for a 100 
GeV DM particle, which for an ordinary DM particle corresponds to a constraint on the DM-nucleon 
cross section of a v SI ~ 6.0 x 10~ 43 cm 2 |93|. This interpretation of the data relies on the assumption 
that DM has come into equilibrium in the Sun, at which point the capture and annihilation rates 
are comparable; see e.g. the review |94|. If the equilibrium is not achieved, the DM annihilation 
rate would be suppressed compared to the capture rate. 

In the DDDM scenario, two factors could enhance the solar capture rate (given the same DM- 
nucleon cross section). First, if the sun is inside the dark disk, DDDM should have a larger local 
density, ~ 10 times as large as the ordinary DM local density near the sun. The capture rate 
could also be enhanced by a larger gravitational Sommerfeld enhancement scaling as the inverse of 
the velocity dispersion |95| , again a factor of 10 compared to that of ordinary DM. The DDDM 
annihilation rate is also Sommerfeld enhanced due to the long-range U(l)i> Thus we expect that 
DDDM always comes into equilibrium in the Sun unless the DM-nucleon cross section is very small. 



The DDDM-nucleon cross section varies from about 10~ 49 cm 2 96 , if DDDM interacts only through 
a loop-suppressed coupling to SM photons, to 10 -44 cm 2 if DDDM couples to SM gluons at the 
one-loop order. This is equivalent to a capture rate in the range 10 17 s^-lO 22 s~\ assuming that 
the DDDM capture rate is 100 times as large as that of ordinary DM with the same cross section 
for scattering on nucleons. Thus DDDM could potentially lead to a signal in the ongoing IceCube 
experiment. It is interesting that, although small velocities make direct detection more difficult, 
they enhance the solar capture rate and could lead to larger signals at IceCube, which then plays 
an important complementary role. 

Though it is not possible in our minimal models, more complicated models could also lead to 
possible signals from high energy gamma rays or charged particles such as near the Sun. For 
instance, DDDM could annihilate into metastable intermediate particles that decay outside the Sun 
into photons or e + ,e~ pairs; or if DDDM particles scatter inelastically, captured DDDM can be 
bound in elliptical orbits of order the size of the Sun and can then annihilate outside of the Sun [97] . 



7.4 LHC Searches 

In principle, the LHC can search for dark matter as well. WIMP searches can proceed if the WIMP 
is part of a larger sector, such as a supersymmetric theory where charged superpartners can decay 
to the LSP. Such searches are unlikely to apply for this new sector unless it is also part of a BSM 
model, which we leave an open question. 



Other searches 98 -100 rely on crossing the interaction (see Appendix |A|) responsible for either 
direct or indirect searches. Ref. |96| studies whether an operator that produces the Fermi signal 
can also lead to a detectable LHC signal. They concluded that the signal is barely detectable when 
there is no large boost factor. Our large boost factor implies a smaller strength matrix element that 
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will not be observable in the near future. 



8 Conclusions and Future Work 

In this paper we have shown that it is possible for a subdominant component with up to 15% of 
all dark matter and dissipative dynamics to collapse into a disk similar to the baryonic one. In 
effect, such matter behaves much like a new kind of ordinary matter, constituting a hidden world 
neighboring our own. If it annihilates to visible-sector particles, we could see a striking enhanced 
indirect detection signal distributed on the sky very differently from that expected for ordinary dark 
matter. Even without such indirect detection channels, it is very possible that such new forms of 
matter could be detected through their gravitational interactions with other matter. We have given 
some simple estimates of the properties of Double-Disk Dark Matter, but much remains to be done. 
Here we will briefly outline some important directions for future work. 

Numerical simulations. Questions such as small-scale structure and the expected alignment 
of the DDDM and baryonic disks would best be answered through numerical simulations of galaxy 
formation (e.g. in mixed A^-body/hydrodynamical codes). A fully correct picture, especially in the 
case that equipartition between heavy and light particles is not fast, would likely require modeling 
dark electric and magnetic fields as well. Such numerical simulations might also shed light on 
the expected velocity distribution of DDDM near the Sun, which is important for understanding 
whether direct detection could ultimately be possible. 

Large-scale structure. DDDM could impact the large-scale structure of the universe in ways 
that might be detectable in the CMB, galaxy/galaxy correlation functions, or other observables. For 
instance, dark acoustic oscillations are a possible signal fl8] . It would be interesting to determine 
whether surveys of large-scale structure, possibly including upcoming 21 cm observations that probe 
the cosmic dark ages, could be sensitive to the existence of DDDM. 

Small-scale structure. ACDM simulations have problems at small scales such as overly large 
cusp predictions and too many satellites. The potential for DDDM to address these problems 
deserves study. After the disk forms and cools, regions within the gas of dark atoms can suffer 
gravitational collapse, perhaps leading to interesting small-scale structure and compact objects. 

Chemistry and nuclear physics. It will also be interesting to study dark matter chemistry, 
which should resemble hydrogen chemistry, as well as nuclear physics if additional interactions are 
included. Further cooling processes that depend on dark atoms and molecules will be interesting to 
study in this case. 

Observational tests for a disk. The Gaia satellite, or other surveys of stars in the Milky 
Way, will study star velocities with unprecedented precision. It is important to see how this can be 
used to map out the distribution of dark matter and test for the presence of structures like a DDDM 
disk through their gravitational effects. Other possible tests could be microlensing from compact 
DDDM objects or lensing from the net effect of the DDDM disk on light from distant objects. 
From the particle physics perspective, models in which DDDM exists but can be only be detected 
gravitationally are conceivable, so it is vital to understand whether current or future observations 
can directly probe its gravitational effects. 

Indirect detection bounds. New analyses of existing data, for instance from Fermi-LAT, 
could be used to set limits on the annihilation rate of DDDM into Standard Model particles. 
These analyses would differ from the standard analysis because DDDM's spatial distribution is very 
different from a typical halo profile. 
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Nonthermal cosmologies. It will be interesting to explore models that generate the dark 
sector asymmetry and possibly relate it to the baryon asymmetry. A more thorough exploration of 
possible scenarios generating a nonthermal symmetric C, C component at late times could also be 
interesting. 

The various issues highlighted above may not all be decoupled. For example, it has been proposed 
that the cusp/core problem is related to supernova explosions that flatten the dark matter cusp into 
a core. This is in tension with the low star-formation efficiency suggested by the missing satellites 
problem [101 . But suppose that the DDDM sector involves violent small scale events that, like 



supernovae, inject energy that could flatten out the cusp. This may not happen in our minimal 
model, but is conceivable in a DDDM model with a closer resemblance to the Standard Model, for 
instance. Such events may not be observable in visible light, and so the tension with the missing 
satellites problem may not exist in this case. This is just one of many possible directions that could 
tie together the physics and astrophysics of DDDM in novel ways. 

We emphasize that, although a dizzying array of particle physics models have been proposed 
for dark matter, most appear from the astrophysical perspective as the same cold, collisionless dark 
matter. Explorations of dark matter with different astrophysical consequences are mostly limited to 
warm dark matter and self-interacting dark matter, with the latter usually assuming a fixed cross 



section for point-like interactions or, at times, velocity dependence 102 -104 ; though see also 105 



Such scenarios with two components of dark matter — the dominant one essentially noninteracting 
and a small component with self interactions — introduce many new possibilities and are as yet only 
weakly constrained. We believe that Partially Interacting Dark Matter and, in particular, Double- 
Disk Dark Matter go far beyond standard scenarios for dark matter and offer very novel prospects 
for dark matter astrophysics. The phenomena we have discussed in this paper could be just the tip 
of the iceberg. 



A Satisfying Kinetic Mixing Constraints in Models of DDDM 
A.l Kinetic Mixing Constraints 



Ordinary photons and dark photons can mix through the operator kF^F^ 74 . Assuming (as 
we have throughout the paper) that dark photons are exactly massless, such a mixing means that 
ordinary matter has a small charge under U(1)d and dark charged matter has a small charge 
under U(1)em- This would keep the DDDM sector and visible sector coupled to dangerously low 
temperatures. In particular, in the early universe, it would allow the Standard Model plasma to 
produce C and C through interactions like e + e~ — > CC or 7 — > CC. When rac ~ 1 MeV this 



will overproduce relativistic C particles and violate the bound from BBN discussed in Sec. 3.1 



unless 106 



k < 10- 9 . (40) 

For even lighter C masses, mc < 10 keV, there is an even stronger bound of k < 10~ 13 to 10~ 14 
from the cooling of red giants and white dwarfs (for details see |106| and references therein). 

Kinetic mixing has the potential to be dangerous because, as a marginal operator, it can be 
generated at any scale. Even GUT-scale particles charged under both groups could lead to violation 
of the experimental bound. Despite this danger there are several ways of avoiding the bound. First, 
we observe that it is consistent to simply set the kinetic mixing to zero, if there are no particles at 
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any scale charged under both U(l)y and XJ(1)d. This is true even if particles charged under the 
two groups interact in other ways, e.g. through exchanging scalar fields. For example, consider the 
possible 3-loop diagram in Fig. [TU] Because both particles and their antiparticles run in this loop, 
the diagram will vanish: the loop on the right with will be canceled by a loop with ip~ . 




Figure 10: A possible 3-loop contribution to kinetic mixing when no particles are charged under both 
sectors, which is zero because it involves both particles and antiparticles. 




Figure 11: Possible 4-loop contributions to kinetic mixing when there are particles charged under both 
U(l)s. At left: this diagram is proportional to Tr(<5 3 Q|,). The condition that TrQ D (Q 3 ) = 0, i.e. that 
the trace of visible charged cubed vanishes in the sector with any given dark charge, is sufficient to make 
this diagram vanish. At right: this diagram is proportional to Tr(Q 2 Qj :i )Tr(QQp). The condition that 
Tiq(Q^ ) ) — is sufficient for it to be zero. 



In fact, kinetic mixing can vanish to high loop order even if there are particles charged under the 
two groups, by choosing charge assignments to satisfy certain anomaly-like conditions. The one-loop 
diagram is proportional to Tt(QQe)), with Q the visible charge of the particle running in the loop 
and Qd the dark charge. If for every choice of Q the sum of the dark charge of particles with 
visible charge Q vanishes, a condition that we denote Tiq(Qd) = 0, the one-loop kinetic mixing 
will vanish. At higher loops there are diagrams with extra photon and dark photon exchanges 
going as Tr(Q m Q' 1 ^) for m,n> 1. An example of such a diagram at four loops is shown at left in 



Figure 11 If we impose the conditions TVq(Qd) = ^QiQo) = ^QoiQ) = ^QniQ ) = 0> an °f 



these diagrams with four or fewer loops will vanish. Another set of diagrams at four loops involves 



two fermion loops, like the right-hand plot in Fig. 11 but these are set to zero by the same trace 



constraints. These are anomaly-like constraints in the sense that they demand certain vanishing 
traces, but they are much more restrictive than anomalies: they apply to scalars in the loop as well 
as to fermions, and they restrict the trace of dark charges in the sector with fixed visible charge and 
vice versa. Although these conditions are restrictive, they could potentially be satisfied, and could 



forbid kinetic mixing up to 5 loops. Such models would be consistent with the bound in Eq. 40 



Another possibility is that U(l)y or U(l)£> is embedded in a nonabelian group. For example, 
suppose that \J(1)d arises from a group SU(2)£) broken by an adjoint Higgs If the lightest 

particle charged under both SU(2)£> and U(l)y has mass M > {&d), we expect that the mixing 



29 



arises from a dark S-parameter operator |107| 



f^^Tr(^^)^. (41) 

This is consistent with the bound k ~ 10 -9 if, for example, SU(2)x> is broken at the weak scale and 
all particles charged under both groups have masses above 10 9 GeV. In such a scenario, the threat 
that GUT- or string-scale physics renders the model inconsistent with the bounds can be avoided. 

Another distinctive scenario is to consider that an unbroken nonabelian dark force remains at 
low energies. In this case kinetic mixing with the photon is completely impossible. One might worry 
that such a force would confine and prevent long-range interactions. However, with the relatively 
small values of cid at which cooling can be effective, the temperature of our dark plasma in galaxies 
will be too high for confinement to occur. Nonabelian sectors that don't confine because they flow 
to infrared fixed points are also a possibility, but in this case we would need more fields with dark 
charge and the BBN bounds on the number of light degrees of freedom would become more severe. 

One final possibility is that the U(l)x) is not exact and the gauge boson is not massless. 

A. 2 Other Interactions Between C, C and the SM Plasma 

In models with indirect detection signals, we assume that the heavy dark sector particles X and X 
can annihilate to Standard Model particles; for instance, we can consider the process XX — > 77. 
If the light particles C have couplings similar to those of X, then the inverse process 77 — > CC 
can produce relativistic C particles at late times. This may be in conflict with the bound on light 



degrees of freedom discussed in Sec. 3.1 Unlike kinetic mixing, this constraint would arise from a 
higher dimension operator, 



1 



CCF^F^, (42) 



which is dimensionally suppressed and imposes weaker constaints than in the previous subsection. 
The cross section is 

<Tt^CC)n^. (43) 

This scattering process potentially keeps C and C in equilibrium with the Standard Model. To 
check this, we compare to the Hubble rate: 

"Wrr -> cc) ggggfig 06q -^T^H m 

H fc~Jl. ~ Uma * A6 • l44j 

V 90#*M P1 

This shows that this process no longer couples C to the Standard Model thermal plasma when the 
temperature drops below 



9 * ~ 250 MeV ( ^— ( — ) • (45) 

0.06Af P J V 200 GeV / VlOO/ V ' 



Thus, if the coupling of X to the dark sector is fixed through mediators with mass of order A at the 
weak scale, 77 — > CC scattering is approximately consistent with the hidden sector decoupling from 
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the visible sector at temperatures above 200 MeV, which we found to be safe in Section [3| Many 
models for the UV completion of this annihilation model are possible, including those illustrated in 



Fig. 12 In fact, for a fermionic model chiral symmetry can be used to further lower the couplings 
of C relative to those of X, rendering the model even safer from constraints from 77 — > CC. 

Notice that for scalar dark matter 0, the operator <jy (pF^uF^ 11 has dimension six and <fi would 
not decouple from the Standard Model until lower temperatures, posing a potential problem for 
AA^ B ^. These models however could be safe as well if the UV completion works as in Fig 13 if the 
coupling to a new gauge boson W' is absent for the C particles (so that they are not produced from 
the Standard Model plasma) but not for the X particles (so that they can annihilate to photons). 



X,C 
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X,C- 



x,c- 



fa 



Figure 12: Possible scenarios for fermionic dark matter giving rise to an annihilation signal in gamma rays. 
At left: s-channel intermediate scalar. (It could also be a spin-one Z' with a final state). At right: box 
topology with charged intermediate states. 



X -< -j----< i'v/wwa 

< x± : 

X > - - - i >- - - - *\z\zKf\z\s- 

Figure 13: Possible scenario for an annihilation signal in gamma rays: scalar loop diagram with a box 
topology. This can give X annihilation and not C production, consistent with bounds, if C does not couple 
to the new gauge boson W'. 
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